Photon Number Statistics of Multimode Parametric Down-Conversion 
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We experimentally analyze the complete photon number statistics of parametric down- conversion 
and ascertain the influence of multimode effects. Our results clearly reveal a difference between 
single mode theoretical description and the measured distributions. Further investigations assure 
the applicability of loss-tolerant photon number reconstruction and prove strict photon number 
correlation between signal and idler modes. 
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The photon statistics of a light beam are the "finger- 
print" of its quantum state, from which a number of use- 
ful measures of non-classicality may be inferred. Coher- 
ent states reflect closest the classical properties of optical 
fields. These states exhibit a Poissonian photon number 
distribution, while sub-Poissonian statistics are a clear 
signature of quantum light. On the other hand, single- 
mode thermal radiation shows super-Poissonian photon 
statistics. However, thermal light sources typically emit 
a multimode field so that their photon statistics are given 
by the convolution of several thermal photon number dis- 
tributions, which, in turn, leads again to a Poissonian dis- 
tribution. In some applications, such as quantum key dis- 
tribution, the photon statistics and their interplay with 
the multimode structure are intrinsically important for 
determining the security of the communication. 

For multi-partite systems another important quantum 
characteristic is the extent to which they are correlated 
beyond what is possible classically. Quantum correla- 
tions, or entanglement, underlie much of the informa- 
tion processing power of quantum systems. For multiple 
light beams it is therefore important to determine the 
correlations between the different modes, and photon- 
detection methods that can reveal these correlations be- 
come paramount. The joint photon statistics of two 
beams can provide a useful measure of correlation and 
the degree of entanglement, also important for applica- 
tions such as decoy state quantum cryptography [l]. 

In this Letter we report, to our knowledge, the first 
direct detection of complete higher order photon number 
correlations and statistics between two light beams that 
are strictly correlated in photon number (a twin-beam), 
using a photon-number-resolved characterization whose 
performance is compromised neither by dark-count noise 
nor by detector inefficiencies. This enables us to detect 
strict photon number correlations for all observed pho- 
ton numbers including the vacuum. We are also able to 
identify features caused by the multimode character of 
the twin-beam source used. 

In general, photodetectors are not intrinsically mode- 
selective, since they have a finite etendue and spectral 
response bandwidth. In this regard they are ideal for 



measurements in the photon- number basis. However, de- 
tectors that are sensitive enough to respond to single pho- 
tons are usually unable to resolve the differences between 
small numbers of incident photons. In contrast, detec- 
tors for large photon fluxes usually have high efficiency, 
but are too noisy to be sensitive to individual photons. 
Hence, there currently exist two main approaches to re- 
cover photon number statistics. 

The first approach, homodyne detection, employs an 
ancillary optical beam, the local oscillator, to increase 
the photon flux level so that it may be registered on a 
standard photodiode. The photon statistics needs to be 
recovered indirectly from a large data set by means of 
quantum state tomography [2 or the use of pattern func- 
tions [3" . Therefore the photon number statistics are only 
recovered indirectly involving the acquisition of a large 
data set. Further, the local oscillator necessarily acts as 
a mode filter, picking out the component of the input 
beam that matches the mode of the ancilla. Hence, the 
multimode character becomes obscured by the detection 
technique. Although multimode homodyne detection is 
possible [4 , it is laborious and requires a very high signal 
to noise ratio to accurately estimate the photon numbers. 

The second approach combines mode-multiplexing 
with single-photon sensitivity to allow for multimode 
state characterization without the need of a local os- 
cillator. However, achieving single-photon sensitivity, 
high quantum efficiency and low dark count noise consti- 
tutes an experimental challenge. The two most common 
types of detectors for such applications are avalanche 
photodiodes (APDs) and visible-light photons counters 
(VLPCs). VLPCs can detect single photons and achieve 
multi-photon resolution by spatial mode-multiplexing. 
They have been used to determine the non-Poissonian 
photon number statistics of degenerate parametric down- 
conversion (PDC), to analyze the statistics of heralded 
number states using twin beams, and to verify the non- 
classicality of the PDC states [5]. While VLPCs exhibit 
very high quantum efficiencies, they suffer from signifi- 
cant dark count contributions and require cryogenic cool- 
ing. Similar difficulties are encountered using sensitive 
CCD cameras at room temperature [6 . Photon number 



resolved measurement can also be achieved combining 
the idea of mode-multiplexing with the ease of opera- 
tion offered by commercial APDs. Here an input pulse 
is distributed in the time domain into multiple temporal 
modes using an optical fiber network £j . In 2006 Achilles 
et al. introduced a new technique for loss-tolerant char- 
acterization of photon number statistics [8 using a time 
multiplexed detector (TMD) to characterize the condi- 
tioned statistics of one- and two-photon states from a 
twin beam. 

Generation of twin beams by means of PDC has be- 
come a widespread way to implement highly correlated 
quantum states. Sub- and super-Poissonian photon num- 
ber statistics [9 , photon anti-bunching ^10^ and quantum 
correlated quadrature amplitudes [11] have been demon- 
strated. More recently, a PDC source was employed to 
prepare "kitten" states by means of photon subtraction 
[12], paving the way for a new approach to conditional 
quantum state synthesis. 

Although some measurements of joint photon number 
statistics of PDC have been performed, a complete and 
accurate measurement of the joint photon number prob- 
ability Pn,m for two beams had not been achieved. Such 
measurements should quantify the photon number cor- 
relations between the signal and idler beams to all or- 
ders, and determine the moments of the unconditional 
(marginal) and conditional photon distributions for each 
beam. Here, we show such a set of measurements. We 
employ two TMDs for detecting simultaneously signal 
and idler beams, for a range of parametric gains. We 
are able to build up the full joint, marginal and condi- 
tional photon number distributions, and to confirm var- 
ious important correlations. The consistency of the re- 
constructed probability distributions across a variety of 
measurements validates our approach, and enables us to 
measure directly the photon number distribution even for 
highly multimode beams. 

Most crucial for loss-tolerant reconstruction of the pho- 
ton number distribution of a twin beam state |^) = 
J^nV^I^'^) ^^ ^ knowledge of losses, or equivalently 
the detection efficiency. Calibration according to the ex- 
tended Klyshko method [3 [8] requires the a priori as- 
sumption of photon number correlation. This assump- 
tion, however, needs verification, and we are able to show 
this experimentally using the TMDs. 

At the heart of our experiment is a TMD, whose 
response [T3] for incoming photon statistics p = 
(po, ' • ' ^Pn^ • ' •) is modeled as 
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The mathematical description includes a loss matrix L(77) 
and a convolution matrix C. The loss relates the individ- 
ual survival probability 77 of one photon to the survival 



of n out of m > n photons. The convolution accounts 
for a probabilistic distribution of n photons into k differ- 
ent time-multiplexed modes whilst taking into account 
an individual population probability P^. Both matri- 
ces act on a photon probability vector p and yield an 
observable click vector g in the detector. The probabil- 
ity matrices L(?7) and C act non-deterministically on the 
measurement outcome on a shot-by-shot basis. Never- 
theless, they provide an analytic and deterministic func- 
tional dependence for state reconstruction when applied 
to an ensemble measurement. For accurate photon num- 
ber measurements, all model parameters 77, P^ and ^need 
to be experimentally accessible. 

Our experimental setup is depicted in Figure [T] A 
diode laser source at 404nm pulsed with a repetition 
rate of IMHz, with spectral and temporal bandwidths 
of 2.0nm and 60ps FWHM respectively, undergoes mode 
cleaning, followed by polarization control. It is sub- 
sequently coupled into a periodically poled KTiOP04 
waveguide chip by means of a focussing graded index lens. 
The nonlinear interaction of pump photons decaying into 
signal and idler photons via a type-II PDC process (at 
central wavelengths A^ = 740nm, A^ = 890nm) results in 
the generation of orthogonally polarized, photon number 
correlated twin beams. 
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FIG. 1: Scheme of experimental setup (see text); measure- 
ment configurations A, B, C and D provide information about 
detector efficiency, marginal and collective statistics, and joint 
correlations respectively. 



After state preparation, a dichroic mirror in combina- 
tion with a red pass filter splits the PDC state from the 
strong pump. A polarizing beam splitter separates signal 
and idler modes, which are then coupled into multimode 
fibers for detection with different APD or TMD measure- 
ment configurations, labeled as A,B,C and D in Figure 
[1] These setups are fiber pigtailed for easy interchange- 
ability with minimal influence on the detection efficiency. 
Tight temporal filtering with precision of 82ps is applied 
during data acquisition for noise reduction. The exper- 
iments are performed at cw-equivalent powers of pump 
light coupled through the waveguide, ranging from 40nW 
to 375nW. 

The different measurement configurations are used for 



a sequence of experiments that allow for checking the con- 
sistency of the recovered information about the state. We 
begin with determining the detection efficiency and then 
use this calibration to analyze the marginal statistics of 
signal and idler mode respectively. The measurement of 
the collective statistics of both modes is then followed by 
the complete measurement of the joint photon number 
statistics of the PDC source. 

Setup A provides the basis for measuring the efficiency 
r]K = ^ from the count rates of coincidence {Re) and 
single (Rs) events. We calibrate the real efficiency r] by 
evaluating r]K at the lowest pump power of 40nW and 
observe signal and idler efficiencies of r]s = 11.7% and 
r]i = 13.7% respectively. Our PDC shows high source 
brightness even at this low pump power, yielding count 
rates of 4950/s and 5764/s for signal and idler. 

In our second configuration, setup B, we study the 
marginal statistics of the signal and idler mode. We 
replace, for example, the idler-APD from the previous 
setup by a TMD in order to obtain the statistics of the 
idler mode. A knowledge of r]i is crucial for the loss toler- 
ant reconstruction of photon statistics. We infer rji from 
a coincidence measurement with the remaining APD in 
the signal mode. We introduce an additional fiber cou- 
pler to run this measurement and find a decline in effi- 
ciency to r]i = 11.3%. We extract all parameters of our 
model directly from the measured data: the TMD click 
statistics ^on a shot-by-shot basis; the efficiency rj from 
coincidences between idler-TMD and signal- APD for loss 
tolerant reconstruction via L; the occupation probability 
Pk of a particular TMD time-mode defining the convolu- 
tion C. The reconstructed photon statistics are obtained 
using a direct matrix inversion of Eq IT] [8 . The resulting 
marginals are shown in Figure [2^ for a range of pump 
powers, displaying several photon statistics with respec- 
tive values Po, • • • 7^4, which correspond to the growth of 
mean photon numbers. Figures of merit reflecting the 
same information are easily computed, such as the prob- 
ability moments (n^), depicted in inset of Figure |2k and 
compared to a Poissonian fit. Reconstructed marginals 
for TMD measurements in signal and idler beams are pre- 
sented in Figure [2)3. The marginals show a good agree- 
ment for all photon numbers between signal and idler, 
which validates our photon number reconstruction as this 
is expected for photon number correlation. The condi- 
tioned statistics (Figure |2]3 inset) are needed for the cal- 
ibration and confirm the preparation of heralded single 
photon states for the lowest pump power as well as the 
good suppression of vacuum contributions for all pump 
powers. 

Using setup C, we extend our testing of the consistency 
of the inversion method. In conjunction with the previous 
measurements, we now compare signal marginals, idler 
marginals and the collective statistics of the twin beams. 
This experiment comprises the collective coupling of both 
signal and idler beams into a single TMD. Thus, the con- 
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FIG. 2: (a) Marginal photon number statistics in idler mode 
for various pump powers; Inset: first probability moments 
(n^) with a fit to a Poissonian distribution, (b) Compari- 
son of signal and idler marginals for lowest and highest pump 
power; Inset: conditioned statistics demonstrate vacuum sup- 
pression, (c) Collective photon pair statistics (see text) for 
lowest and highest pump powers. 



figuration renders a strict test for the occurrence of pho- 
ton pairs. The suppression of odd photon numbers, pre- 
sented in Figure [2J3, provides direct evidence for photon 
number correlations. Clearly, this configuration cannot 
be used to measure coincidences between signal and idler 
photons. We therefore used the calibration of r]s and r]i 
from previous measurement configurations. As expected, 
the photon statistics exhibit oscillations in the collective 
probability distribution, arising from the pairwise gen- 
eration of photons. Further, comparing the collective 
statistics to the marginal distributions of Figure (2)3 we 
can prove that the distributions are consistent across all 
measurements, because P2n^^^*^^^ ~ p^^^^^ = p^^®^ holds. 
We complement our state characterization with 
setup D and demonstrate the first direct joint photon 
number measurement by simultaneously operating two 
TMDs, one for each beam. This enables us to measure 
directly all photon number correlations between signal 
and idler. The theoretical model of a single TMD can 
easily be extended to cover multiple TMDs. Let us as- 
sume that signal and idler modes obey a photon number 
distribution for which Pn,m is the probability to find n 
and m photons in signal and idler beams respectively. 
Each beam is analyzed by one TMD independently: 



Qk,l = {Cs'LsiVs) ^ CiLi{r]i))pn,r. 



(3) 



Experimentally we find that the probability correlation 
matrix Pn,m is dominantly populated on the diagonal. 



which proves the strict photon number correlation be- 
tween signal and idler. This, in turn, is a strong jus- 
tification for exploiting the a priori assumption of pair 
generation for our source, used for calibrating the detec- 
tor efficiency. In this experiment we coupled the signal 
beam into a single mode (SM)-fiber TMD, and the idler 
beam into multi mode (MM) -fiber TMD. The electronics 
records correlated click statistics between both TMDs. 
The coupling into SM-fibers comes along with a decrease 
of detection efficiency to a level of 77^ = 2.74% whereas 
rji = 11.1% remains unchanged. Our setup is sufficiently 
loss tolerant, though, to cope with the low efficiency. 
The result Pn,m for a measurement with pump power 
of 220nW is depicted in Figure [3| Slight uncertainties in 
the measured click statistics due to finite acquisition time 
of 300s and the modest efficiency in the signal mode con- 
tribute to minute errors found in off-diagonal elements. 
From this data we determine a normalized covariance of 
0.996 in the photon number statistics, compared to 0.052 
in the raw data. The reconstructed twin beam photon 
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FIG. 3: (a) Joint photon statistics Pn,m measured with two 
TMDs and pump power of 220nW. (b) Marginal statistics of 
idler with pump power of 375nW. The measurement is fitted 
against thermal and Poissonian distributions. 



At high power levels we find a mismatch between 
the observed photon statistics and the predictions from 
a single-mode theory of the PDC process, equally de- 
scribed as a two- mode squeezer. The evaluation of the 
quadratic Hamiltonian gives for the signal and idler mode 



a thermal marginal distribution Pn 



(nV 



In con- 



(l+(n))-+i' 

trast, our results obey a Poissonian marginal distribution 
Pn = e~^^^^^i-, as depicted in Figure 3 and illustrated 
by the fits. The thermal fit misses the observation in 
the one-photon contribution by 0.0614, which is more 
than seven times larger than the statistical error. This 
error is derived analytically by Gaussian error propaga- 
tion applied to Eq. [T] taking into account uncertainties 
in each model parameter. The measurement uncertain- 



ties for the efficiencies are 0.009, and for the click statis- 
tics a conventional square root dependence is assumed. 
The discrepancy between the observed Poissonian statis- 
tics and the thermal distribution of a two-mode squeezed 
state can only be explained by a PDC process, introduc- 
ing a multitude of different modes. The convolution of 
many independent thermal statistics yields the observed 
Poissonian statistics. 

In summary, we have demonstrated the first explicit 
measurement of correlations and complete joint photon 
number statistics of a PDC source. Our approach to use 
a TMD for loss-tolerant state characterization enabled us 
to study the marginal photon number distributions in de- 
tail. This revealed a clear signature of a multimode struc- 
ture of the source. We expect the multimode structure of 
light to play an increasing role for the development of ad- 
vanced quantum key distribution protocols as it could be 
exploited for implementing higher dimensional informa- 
tion coding [14 . Multimode states need to be considered 
for developing hybrid quantum communication schemes 
that include continuous variable states as well as single 
photons [15]. 
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